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Abstract 

Finite XXZ chain with double boundaries is considered at critical regime — 1 < 
A < 1. We construct the eigenvectors of finite Hamiltonian by means of the vertex 
opertors and the quasi-boundary states. Using the free field realizations of the 
vertex operators and the quasi-boundary states, integral representations for the 
correlation functions are derived. 



1 Introduction 

The XXZ chain is a fundamental model in understanding of the integrable systems. 
Many attentions have been paid to the XXZ integrable systems [l], [|. The purpose of 
this paper is to derive correlation functions for finite XXZ chain with double boundaries 
at critical regime — 1 < A < 1, by means of the free field approach. 

In the earlier works || |4[] the XXZ chain with a boundary was considered at massive 
regime A < —1, in the framework of the free field approach. The integral representa- 
tions of the correlation functions were derived. It was shown that boundary quantum 
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Knizhnik-Zamolodchikov equations with the certain shift, governed the correlation func- 
tions. The f/ g (s/ n )-generalization of the papers on the XXZ chain [Q, |J was given in 
||, ^j. In the paper results for finite XXZ chain at massive regime A < — 1 were 
extended to critical regime — 1 < A < 1, using bosonizations of vertex operators M. 

Y.Fujii and M.Wadati |J noticed that solutions of boundary quantm Knizhnik- 
Zamolodchikov equations without shift became eigenstates of finite XXZ chain with 
double boundaries. They constructed eigenstates of finite XXZ chain with double 
boundaries at massive regime A < —1, by means of the free field approach. In the 
present paper we shall study finite XXZ chain with double boundaries at critical regime 
— 1 < A < 1. We shall derive the correlation functions as integrals of meromorphic 
functions involving Multi-Gamma functions. 

Now a few words about organiozation of this paper. In section 2 we formulate the 
problem. In section 3 we construct the realizations of eigenstates. In section 4 we derive 
integral representations for the correlation functions. In Appendix A we summarized 
the bosonizations of the vertex operators ||. In Appendix B we summarized the Multi- 
Gamma functions. 



now 



2 Boundary quantum KZ-equation 

In 1984 I.V.Cherdnik [10] proposed the following systems of difference equations, 
called boundary quantum Knizhnik-Zamolodchikov equations. 

F{Pi, ■ ■ • ,&+iA,-- - ,p N ) 
= T i (A,---,^|A)F(/3 1 ,-..,/3 i) -.. ) /3 jV ), (j = l,-..,N), (2.1) 

where the shift operator Tj(@i ■ ■ ■ /3jv|A) is given by 

Tj(Pi, ■ ■ ■ ,/3at|A) = R jd _ 1 (J3 j -/3 j - 1 + i\)---R jtl (J3 j -/3 1 + i\)K j (J3 j ) 

x Ri t j(Pi + (3j) ■ ■ ■ Rj-i t j(Pj-i + (3j)Rj +1 j(Pj + i + f3j) ■ ■ ■ Rn,j(Pn + Pj) 
x K^R^iPj-p^.-.R^ifr-p^). (2.2) 

The i?-matrix R(P) and the boundary i^-matrix K(P),K(P) are specified later. The 
solutions of the boundary quantum KZ equations represent various physical quantities. 
For the case of the shift parameter \ = 2n, the certain solutions of the quantum KZ 
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equations represents iV-point correlation functions for the massless XX Z chain with a 
boundary [|7|, which is described by the following Hamiltonian : 

1 oo 
n=l 

where we set a parameter — 1 < A < 1 and a parameter h represents the boundary 
external field. The a®, and cr* stand for the Pauli matrices acting on the n-th site of 
the half Infinite spin chain : ■ ■ ■ ® C 2 <S> C 2 ® C 2 . The author [^] derived the integrate 
representations of the correlation functions for the above model. 

In the present paper we shall consider the case of the shift parameter A = 0. In this 
case the solution of the boundary quantum KZ equation fl2.1| ) represents the eigenvector 
of finite XX Z chain with double boundaries at critical regime — 1 < A < 1. The 
Hamiltonian Tip of our considering model is given by 

1 N ~ l 

n F = --J2 {<+i< + °l+x< + + h x a{ + h N a z N , (2.4) 

n=l 

where we set a parameter — 1 < A < 1. Parameters hi,h^ represent the boundary 
external fields. The a®, and a z n stand for the Pauli matrices acting on the n-th site of 
the Finite spin chain : (C 2 )® N . 



Let us set the i?-matrix as 



/ 



R{p) = r{p) 



\ 



b{(3) c{(3) 
c(P) b(P) 



(2.5) 



where we set the components as 



sh 







sh 



/3 + iri 



Here we set 



sh 



711 



sh 



S 2 (lp\27l, 7T(e + l))S 2 (-i/3 + 7T|27T, tt{£ + 1)) 

'S 2 {-iP\27T, n(t + l))S 2 {%0 + 7r|27r, tt(£ + 1)) ' 



(2.6) 



(2.7) 



3 



where S 2 {P\ujiuj 2 ) is the double sine function defined in Appendix B. 

Let {v + i V-} denote the natural basis of V = C 2 . When viewed as an operator onV ®V, 

the matrix elements of R(f3) are defined by 



h J2=± 

The i?-matrix satisfies the Yang-Baxter equation : 



(2.8) 



Rn(Pl ~ fo)Rl3(Pl ~ Ps)R2 3 (P2 - Ps) = R23(P2 - Ps)Ri 3 (Pi - Ps) R^Pi - P2) ■ (2.9) 
The normalization factor tq{0) is so chosen that the unitarity and crossing relations are 



Rn(P)R2i(-P) = id, 



R{-(3) 



kik 2 



(2.10) 
(2.11) 



Let us set the boundary ^-matrix K((3) by 



K(0) = k(J3) 



sh 







\ 

v + P 
^ + 1 



V 



sh 



v-P 



where the normalization factor is given by 



(2.12) 



k(P) = ko(P)h(P), 



(2.13) 



where 



UP) 



S 2 (-2iP + 47r|4vr, vr(£ + l))S 2 (2ip + 37r|47r, tt(£ + 1)) 



HP) 



S 2 (2ip + 4vr|4vr, vr(£ + l))S 2 (-2iP + 37r|47r, tt(£ + 1)) ' 
S 2 {-ip + + 7r\2n, tt(£ + l))S 2 {ip + + 2n\2n, n{£ + 1)) 
S 2 {ip + iu + 7r|27r, tt(£ + l))S 2 (-ip + iu + 27r|2vr, vr(£ + 1)) ' 



The matrix elements K(P) k are defined by 



(2.14) 
(2.15) 



K{P)v k = Y. v i K ^) k 3- 

j=± 

The i?-matrix and the if-matrix satisfy the Boundary Yang-Baxter equation. 



(2.16) 



K 2 {p 2 )R 2l {p x + p 2 )K 1 {p 1 )R 12 {p 1 - p 2 ) = R 2l {p x - P^MR^ + p 2 )K 2 (p 2 ){2.17) 
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The normalization factor k((3) is so chosen that the boundary unitarity and the boundary 
crossing relations are 



K((3)K(-(3) = id, (2.18) 

K (p + f)[ = g^)^(-^ + y)^ ( 2 - 19 ) 

Let us set the boundary i^-matrix K(0) by 

K(f3) = K(f3)\^ u . (2.20) 

Note. For the another shift parameter A = 2tt case, we take another choice of the 
K -matrix K(j3). See the reference ffij. 

The derivatives of .R-matrix and i^-matrix are given by 



—R jJ+1 (P)P jjj+ i 



d_ 

2^1) x Ad \ (r '^'- + ^ + A(T ^ +i) + consL (2 - 21) 



sh 



d_ 

d(3 



Km 



ch 



1 „ u + i 



x >^ — (-a* + const. (2.22) 



/3=0 S ■ gh 



Here the anisotropic parameter A = —cos 

Therefore the shift operator Tj(f3\ ■ ■ ■ /3n\0) is related to the Hamiltonian Tip (2.4) as 
following. 

sh (jTl) x x ( Jr Ti ) (0 ' ' ' ' ' 0|0) = Hf + consl (2 ' 23) 

Here the boundary magnetic fields hi, are related with parameters /i, v, 



^ = -i x ah ( ) + ^ = -ixBh(-g-) ^ + (2.24) 



e+iy ve + i 
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We have 

Tj(0,--- ,0|0) = id. (2.25) 
Let us set the eigenvector \/3i ■ ■ ■ (3^) by 

Tj(Pi,--- , P N \0)\Pi, ■ ■ ■ ,p N ) = ,Pn). (2-26) 

Let us set the dual eigenvector by 

{fix, ■ ■ ■ ,Pn\Tj(Pi, ■ ■ ■ ,Pn\0) = (Pi, ■ ■ • ,Pn\- (2-27) 
The above eigenvectors satisfy the followings. 

H F \0,--- ,0) = Const.\0,--- ,0), (0, • • ■ ,0\H F = Const. (0, ■■ ■ ,0|. (2.28) 
In the next section we shall construct the eigenvector \/3i ■ ■ ■ Pn) explicitly. 

3 Eigenvectors 

In this section we solve following eigenvector problem. 

%{Px,--- , Pn\0)\Pi, ■ ■ ■ ,Pn) = \Pi,---,Pn). (3.1) 
The eigenvector is realized by using the vertex operators $-,(/?). 

\Pu---,Pn) = T^C J2 (G\$eM---®e N (PN)\G)(v ei ®---®V eN ). (3.2) 

ei---ejv=± 

Here the vector \G) and the dual vector (G\ are characterized by the following relations. 



(G|%(/3) = K{P)i{G\^{-p), (j = ±), (3.3) 
^(-P)\G) = K(P)^(P)\G), (j = ±), (3.4) 

We call the auxiliary states (G\ and \G) "quasi-boundary state". Using the following 
commutation relation of the vertex operator and the characterizing relations (|3.3|. p.4|) 



of the quasi-boundary state, we have the equation (|37T 



Q jl (p 1 )*M= E R(Pi-P2) k j 1 J:$kM®M- (3-5) 

ki,k'2=± 
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We give the free field realization of the quasi-boundary state. 
Let us introduce the free bosons b(t), (t G IR) by 



sh I — ) sh(7rt)sh 



tsh ( — 

V 2 

Let us set the Fock space 7i generated by the vacuum vector (vac\ which satisfies 

(vac\b(-t) = 0, if t > 0. (3.7) 
The quasi-boundary state (G\ is realized as followings. 

(G\ = (vac\e G , (3.8) 

Here we have set 

G -2j [b(t),b(-t)] m dt + J Ut),b(-t)] bm ' ^ 

where 

G 2 (t\fx) = -1, (3.10) 
1 sh C^t) sh (H/i + |flf) i sh (jt) sh (^) ch 

Let us prove the relation (|3.3|). In what follows we use the abberiviations. 



= exp (- J™ ^/ Pt d^j , U-(0) = exp QT ^e^dt^j , (3-12) 



oo 



U+{a) = exp ( / j^e^df I , I/L(a) = exp I - / e — ) . ( 3 .13) 
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In what follows we omit non-essential constant factors. 
At first we explain the formulas of the form 

X(A)F(/5 2 ) = Cxyifa - fa) : X{f3 1 )X{(3 2 ) :, (3.14) 

where X,Y = Uj, and Cxy{@) is a meromorphic function on C. These formulae follow 
from the commutation relation of the free bosons. When we compute the contraction of 
the basic operators, we often encounter an integral 



OC 



F(t)dt, (3.15) 



o 



which is divergent at t = 0. Here we adopt the following prescription for regularization 
: it should be understood as the countour integral, 

where the countour C is given by 



Contour C 

The actions of the basic oprtaors on quasi-boundary state (G\ are evaluated as followings. 

(G\U-(J3) = Const.m{/3)(G\U + {-(3), (3.17) 
(G|E/_(a) = Const. J(a)(G\U + (-a). (3.18) 



Here we have set 

m(j3) 



T 2 (2if3 + 47r|2a^2)T2(^ + tt(£ + l)]2a^) 
T 2 (2z/3 + 37r|2a;ia;2)r 2 (22/? + tt(^ + 2)|2wicu 2 ) 
T 2 (i(3 + ifx + 7r\uiuj2)T2(if3 — + 7r(£ + 2)|cjicj 2 ) 
r 2 (i/3 + z/i + 27r|u;iu; 2 )r 2 (i/3 - i/i + 7r(£ + l)^^) 



we have proved the "+"-part of the characterizing relation (!T3) 



(3.19) 



■p / — i)i+ia _|_ -1 __ 1 
V 7r(f+l) ~ l ~ 2(£+l) , 

J(a) = a x v / ^iL. (3.20) 

■p / ifi+ia _|_ 1 \ 

We have 

(G|$ + (/3) = m(0)(G\U+{P)U + {-0). (3.21) 
Because the function ra(/?) satisfies 

*<flt = ^ (3.22) 



We will prove the "—"-part of the equation ( |3.3| ). Using the actions formulae of the basic 
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operators on the quasi-boundary state, we have the following. 

(G\*-(J3) = Const.m(P) fj aXaX U T + if } + 2fcTT)) 



ei,£2=± 

a + P ni \ 1 I 7r(£+i) ~t 1 2(5+1) 



x sh - — - + 



£+1 2(f+lW F f i/H-ia i 1 

x (G\U + (J3)U + (-p)U+(a)U+(-a). (3.23) 

Note that the operator part of the above equation : (G|[/_(/9)C/_(— (3)U_{a)U-(— a) is 
invariant under the changes of variables : a <-> — a,/3 «-> — /3. 
We have 

m^sh - m(-(3)-hh (|±|) (G|$_(-/3) 
= Const, x sh — - — / da T + — 



«(— /U + eat) 1 \ T-r /yU + ea m 



TT r , / + 1 -— — - x a TT sh 



_ : x tt(£+1) 2(^ + 1); 1A V € + 1 2(^+1) 

x (G|?7_(/3)C/_(-/3)0"_(a)C7_(-a). (3.24) 

The integrand of (RHS) is anti-symmetric to a change of integral variable a «-> —a. It 
means the left-hand side becomes zero after taking integral. Therefore we get 



m 



-/3)sh (l^) (G|$_(/3) = m(/3)sh (|±|) (C?|$-(-/3). (3.25) 



We have proved the "—"-part of the characterizing relation 
The quasi-boundary state \G) is given by the following. 



\G) = e G >ac), (3.26) 



where 



G*- 1 r b( t?dt+ r °* m b( t)dt (3 27) 

G -2j MA=W} { )d+ Jo W),K-t)r > ' { ] 



where 



G*(t\») = G 2 {t\/i), GKtl/i) = -G^u). (3-28) 
As the same manner as the above we can prove the characterizing equations (|3.4j) . 
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Note. When we consider the massless XXZ chain with a boundary jffl. We introduce 
the boundary state \B) and the dual boundary state (B\. Left quasi-boundary state (G\ 
differes from the dual boundary state (B\. Right quasi-boundary state \G) coincides with 
the boundary state \B). Physically the boundary state \B) of the paper ^ coresponds to 
the vacuum expectation value ■ ■ • $>((3n)\G) of the present paper. Both quantities 

\B) and (G\§((3i) ■ ■ • $(/?at)|G) represent an eigenvector of the Hamiltonian for each 
model. 

Let us construct the dual stationary state, 

Av|2}(A,-- - ,Pn\0) = (Pi,---Pn\. (3.29) 
The dual stationary state is realized by using the dual vertex operators 5>^(/?). 

(/v--Avi = -^y £ ( f \K(pj---kam\f)K®---®< n )- (3.30) 

ei,— ,ejv=± 

The dual vertex operators are related to the vertex operators, 

= <M/3 + iri), (j = ±). (3.31) 

The quasi-boundary state (F\ and \F) are characterized by 

(F\^(i3) = KffljiFWi-P), (j = ±), (3.32) 
= K((3y^*((3)\F} } (j = ±). (3.33) 

The quasi-boundary states (F\ and \F) are realized as followings. 

(F\ = (vac\e F , \F) = e F *\vac). (3.34) 

Here we have set 



2 Jo W),K-t)] Jo [K«),K-«)1 

- 2I mM=m b{ ~ t)dt+ l wfM-M (3 ' 36) 

where 

F 2 (t\n) = -e- 27rf , (3.37) 



-(y)-(?) ch (K 



sh(|(e + l)t) * sh + 



(3.38) 
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F*(t\u) = -e 2rf , (3.39) 
F*(t\u) = -e 2nt xF^tlu). (3.40) 

The characterizing relations are proved as the same manner. Here we omit details. 
Formally other eigenvectors are constructed by inserting the type-II vertex operators. 

E ( G \*M) ■ ■ •^(Mntfi) • • • ^Mm)\G)(v 61 ® • • • ® v eN ). (3.41) 

ei— ejv=± 

4 Correlation Functions 

In this section we calculate the vacuum expectation values of the type-I vertex oper- 
ators, and obtain them as integrals of meromorphic functions involving Multi-Gamma 
functions. We compute the following 2/V-point function, 

P^.., eW ,.., ei;?? ({/3*}; {&■}) 

<^l^l(ffl-"^(^)l^> x (G v \^ N ((3 N )--^ €1 (f3 1 )\G v ) ^ (41) 

Here we set the state (G v \, \G V ) and (F v \, \F V ) by 

(G v \ = {vac\e°\ \G V ) = e°^\vac), (4.2) 



(F„| = (vac\e F \ \F V ) = e F ^\vac). (4.3) 



where 



G , = t^m Ktfdt+ rjm-w, (4.4) 



2j„ [6(t),K-«)l Vo [K0.H-01 

i r° e -"G;(t| M )) 2 r gjjfjrt 

= U [KO.K-01 i <-')* + y M^MI 6 *-')* (4 ' 5) 

f v = 1 r^M 4(t) . (ft+ rjm_ bm , ,4.6) 



2i [6(0,6(-0] v ' Jo IOO.M-01 

We have 

WGJ = <G|, hm|G,) = |G), limjFj = <F|, lim|F,) = |F). (4.8) 

First we compute the normal-ordering of the vertex operators. 
Let us denote by A the index set 

A = {a\e a = -,1 < a < N} (4.9) 
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We have the following expressions. 

(G ri \^ N (P N )---^ 1 (Pl)\Gr l ) 



(G v \G ri ) 

r 2 (t(p b2 -p bl ) + 2tt|ccw 2 ) r 2 (i(p b2 -p bl ) + ^ + i)Kn>2) 

, , ; v T 2 (i(p b2 - (3 bl ) + 7r|wiw 2 ) r 2 (i(A, 2 - /3 6l ) + tt(£ + 2)|^w 2 ) 



n 



oo 



i(«a - Pa) , 1 \ ^ / -«„) 1 



x da a [[T( + ___)r[ + 



7r(e+i) 2(e + i); v 2($+i) 



X 



tt I , v V 7r(£+l) £+1 

n, „„C 4 I 1 I / ^ . -, \ 




X 



»i.« 2 eA ^ V ^(^+1) £+1 

, , ^(A> - go) 1 \ p ( Ma - Pb) 

1 7r(e+i) 2^+1); tt v 1 2«+i) 



tt v 2: ^ ^ i til TT 



aS A,i< 6 <iv ^tt^+i) 2 (f + I) J ^ A ^ N \tt(Z + 1) 2(^+1) 

x I v ({(3 b }\{a a }). (4.10) 

Let us denote by A* the undex set 

A* = {a\e* a = +,1 < a < N} (4.11) 

We have the following expressions. 



n 



r 2 (i(A 2 -p bl ) + gT^wg) r 2 (i(& 2 - /? 6l ) + 7r(e + lj^wg) 

, ; v r 2 (i(/3 62 -/3 6l ) + 7r|c; 1 c; 2 )r 2 (i(/3 62 -/3 6i ) +7r (£ + 2)1^2) 



X 



ax<a 2 
a l i a 2 ^ ^* 



x 



n 




r , i(Pb -a a ) 1 \ r / i(Q« - Pb) 1 



7r(£ + l) 2(^+1)7 TT W£ + l) 2(^+1) 



a<b A (A - a,) 1 1 \ / j( ao _^) 1 

" eAW ^((+1) + 2(e + 1); aeA *' 1 ^ JV V tt(C + 1) 2(^ + 1) 

x ttOIK})- (4.12) 
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Here we have set 



T(SRXU u (ggl exp CC exp {£ Y A (t)b(t)dt) \G V ) 

W - - W <r XMm - t: ^ ^ ^ YMdt) |F -> . (4.14) 

\ F v\ F v) 



with 



X A (t) = y4^-y^TT, (4-15) 



X gift,* gictai 



F A (t) = + (4.16) 

We evaluate the quantities I v ({/3b}\{&a}) , In{{@t ~ ^}|{«a})- Using the completeness 
relation of the coherent state [§, [7j , and performing the integral calculations, we have 

W}|{«.}) 

/ r° 1 sh (f ) sh (Trt) sh (?f ) 
eXP ^7 1 - e- 2 ^ tsh (f(f + l)t) 

x ^-i e -^X A (i) 2 + e- 2 " t X A (t)F 4 (t) - ie-*y A (t) a ) (4.17) 
' {(^(tl/x) - e^Gt(t|z/))X A (t) + (GI(t|i/) - e^G^/i))!^)}^ 



o 1 - e~ 21 ?* 
and 

-ttOik}) 

= exp 



1 sh(f)sh(7Tt)sh(^) 



1-e" 2 "* £sh(f(£+l)t) 
x ^-I e -^- 2 -*X^(t) 2 + e- 2 *X A *{t)Y A *{t) - ^ t+27Tt Y A *(t) 2 ^j (4.18 

+ J°° ^T-^imiri - e-^F*(tk))X A ,(t) + (F*(t\u) - e^^F 1 (t\^))Y A .(t)}d^ 
In what follows we use the abberiviations : 

uji = 2ir, lu 2 = 7r(£ + 1),lu 3 = 2rj, fi + = //, /i_ = v, (4.19) 
The vacuum expectation value is evaluated as following. 
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Here we set 



TTTT / S 3 (2ief3 b + it + rj\uiuj 2 uJz) tt tt r 2 (ae(/3 bl - (3 b2 ) + 7r|^i^ 2 ) 
ii H V S ^Pb + 27r + ^1^3) 11 r 2 (ie(p bl - (3 b2 ) + 2 7 r|o; 1 a; 2 ) 

TT TT g3(H^bi + + 7T + g |^1^2^ 3 ) ^3 ~ A»a) + TtIo;^^) 

6 <1 i± ^(^(/^ + A, 2 ) + 2vr + TjluxWaWa^ieOflfe - A> 2 ) + 2n\u 1 u 2 u 3 ) 

tt -TT r 3 (ie/3 fc + ijXe + Trla^a;^^^/^ - + n£ + 27r|c^ 1 c^ 2 ^3) 
Ai- -LA r 3 (ie/5 b - i/i £ + 7r£ + 7r|a;ia; 2 a;3)r3(ie/3 6 + i// e + 27r|u;iu; 2 u; 3 ) 

tt jr r 3 (— ie/?6 + V + ifi e + 7r|c^i^2^ 3 )r 3 (— ie/?b + i] - i[i e + n£ + 271^0;^) 
AA £ A1 Y 3 (-ie(3 b + r] - ifi e + ir£ + 7r\uiU! 2 uj3)T 3 (-ie(3 b + r] + i/j, t + 27r|a;ia; 2 a;3) 

N 



nn 

6=1 e=± 
N 

nn 

6=1 e=± 
N 



Y 3 (2ief3 b + 7r|2cui, cj 2 , 2w 3 )r 3 (2ie/^ + 47r|2a;i, u; 2 , 2cu 3 ) 
r 3 (2ie/3 fe + 27r|2a;i, cu 2 , 2w 3 )r 3 (2ie/5 b + 3^2^, w 2 , 2a; 3 ) 

r 3 (2ie/3 6 + tt£ + tt|2wi, w 2 , 2cj 3 )r 3 (2ie/? 6 + ?r£ + 4tt|2wi, w 2 , 2w 3 ) 



r 3 (2ie/3 6 + < + 2tt|2wi, cj 2 , 2cj 3 )r 3 (2ie/5 b + vr£ + 3tt|2wi, cj 2 , 2cj 3 ) 
r 3 (2ie/3 6 + 2?7 + 7r|2wi, u; 2 , 2u 3 )T 3 (2iep b + 2r] + 4n\2w u w 2 , 2^ 3 ) 



f = i M V r 3( 2?e ^ + 2r ? + 27r l 2 ^i' ^ 2 ,2w 3 )r 3 (2ie/3 6 + 2 V + 37r|2^, cu 2 , 2w 3 ) 1 ' J 

yr l r / T 3 (2ie(3 b + 2r] + ir£ + 7r|2a;i, o^ 2 , 2cj 3 )r 3 (2ie/j b + 2?y + jrj + 4^2^, u 2 , 2uj 3 ) 
f = i M V r 3( 2ie ^ + 2r ? + < + 2tt|2wi, cu 2 , 2cu 3 )r 3 (2ie/5 fe + 2 V + vr£ + 3tt|2wi, cj 2 , 2w 3 ) 



x 



X 



X 



OK}) 




,S'o(2ieo! a + r/|^ 2 ^ 3 ) -q J~J \ " (4+u ^ T 



i(a ai -cta 2 ) 1 



S 2 (2iea a + 7r£ + t]\uj 2 uj 3 ) 11 A ^ r f 'Kr^) , i N 

ai<a 2 e=± 1 ^ tt(£+1) ^ , 

g 2 (^(Q!ai + «a 2 ) + r?|a; 2 cc; 3 )5 , 2(i6(a ai - a a2 ) + tt|^ 2 u; 3 ) 
± S 2 (ie(a ai + a a2 ) + 7r{ + 77|a; 2 cc; 3 )S , 2 (ie(a ai - a a2 ) + ?r(£ + l)!^^) 

II II - 2 ^ aa + ^ + ^ _ i ^ UJ2UJ3 "> T2 ^~ i(Eaa + ^ + ^ + f - ?/ie|^2^3) 



nn 



a£A e=± 



r 2 (iea a + f + iyu e |a; 2 a; 3 )r 2 (-iea a + r] + | + i/z e |cc; 2 ci; 3 ) 




r 2 (2iea a + 7r 


^ 2 ,2cu 3 )r 2 (2iea a + 7r(^ + 1) 




r 2 (2iea a 


a; 2 , 2uj 3 )T 2 (2iea a + tt!;\uj 2 , 2uj 3 ) 



aeA e=± 



r 2 (2tea a + 2r] + tt|cj 2 , 2a; 3 )r 2 (2tea a + 2?? + ?r(£ + l)\cu 2 , 2u 3 ) 
T 2 (2iea a + 2r]\u 2 , 2u 3 )T 2 (2iea a + 2rj + 7r^|cu 2 , 2a; 3 ) 
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#*(WIK}) 



TV 



S 2 (i(a a + (3 b ) + 77 + vr£ + l\u 2 u 3 ) 



TT TT ^IH^a TWT')T/n,T g l^^ 3 ; , s 



5 2 (i(o; a + /3 6 ) +77 + f|a;2^3) 

nnnM%^+^H(^-A)+* 



X 

aeA 6=1 e=± 



^2^3 



The vacuum expectation value is evaluated as following. 

i;m - «}iK}) = ifmw a m}\{* a })i; a ({a a }). (4.24) 

Here we set 



1 r l r / Steffi + 7T + 27re + 77^1^3) tt tt r 2 (ie(/3 6 * - + 7r|wiW 2 ) 

1111 \/ RJOieR* _|_ Ott + 9-71-^ + nLiw.w,^ H H 



- ± V 5 3(2ie/3 6 * + 2tt + 27re + r^u^) ^ e A = A r 2 (ie(/3 6 * - /%) + 2tt|^ 2 ) 
tt tt S 3 (te(P* bl + /%) + tt + 2tt6 + t?)^^)^^ ~ /%) + ^1^3) 
b J 2 e =i ^(^(/^ + /%) + 2tt + 27re + 7 7 |u; 1 u; 2 u; 3 )S 3 (7;e(/3* i - /%) + 27r|^c^ 3 ) 

tttt r 3 (ie/3£ - «A*e + vr^ + 27r + TrelwiU^u^jTs^e/^ + i// e + 7re + nlou^Us) 
11 11 r 3 (ie/3 6 * + i/x e + 7re|cc;iCt;2u; 3 )r 3 (ie/3^ - i/i e + tt£ + 3tt + ire^u^) 

TV 

X 

6=1 e=± 
TV 

X 

6=1 e=± 
TV 

X 

6=1 e=± 

iV 

X 

6=1 e=± 

TV 

X 

6=1 e=± 
TV 

X 

6=1 e=± 



nn 

6=1 e=± 

TV 

nn 

6=1 e=± 

TV 

nn 

6=1 e=± 

TV 

nn 

6=1 e=± 

TV 

nn 



r 3 (— ie/3£ - i// e + 77 + 7r£ + 27T - 7re|u;iu;2u; 3 )r 3 (— ie/^ + 77 + i/j, e + tt - Tie\uJiUJ 2 uo 3 ) 
T z (-itPl + + 77 - 7re|cJiCJ 2 cj 3 )r 3 (-7e / 95 - i/i £ + 77 + ir£ + n - 7ie\uJiUJ 2 uj 3 ) 

T 3 (2ie(3; + 2ire + tt^x, u 2 , 2u 3 )T 3 (2ie(3; + 2ire + 47r|2o^ o; 2 , 2u 3 ) 
r 3 (2ie/3* + 27re + 2^2^ u; 2 , 2w 3 )r 3 (2ie/? 6 * + 27re + 3^2u^ u 2 , 2u 3 ) 

T 3 (2ie(3* + 2ire + tt£ + jjq^a^, 2o; 3 )r 3 (2te^ + 2?re + ?r£ + 4tt|2^i, cu 2 , 2cj 3 ) 
r 3 (2ie/3* + 2ne + vr£ + 2ir\2u u u 2 , 2u 3 )T 3 (2iep; + 2ne + tt£ + 3ir\2u u u 2 , 2uj 3 ) 

T 3 (2ie(3* + 2ire + 2r] + tt|2wi, oj 2 , 2uj 3 )T 3 (2ie(3* + 2ire + 2r] + 4tt|2wi, co 2 , 2lo 3 ) 



T 3 (2ie(3* + 27ie + 2 V + u 2 , 2u 3 )T 3 {2ie(3* + 27re + 2 V + 3vr|2cui, u 2 , 2lj 3 ) 



V 3 {2ie(3* h + 2ne + 2r] + tt£ + %\2tu u u 2 , 2uj 3 ) 
r 3 (2ie/3* + 2ne + 2rj + tt^ + 2ir\2u u cu 2 , 2co 3 ) 



TTTT / T 3 (276^ + 2716 + 2 V + < + 47r|2^, UJ 2 , 2uj 3 ) 
1111 V r 3 (2ie^ + 27re + 277 + 7r^ + 37r|2wi,w 2 ,2w3)' 
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N 



X 



nn 

aGA* 6=1 
N 

nnnu 

a£A* 6=1 e=± 



S 2 (i(a a + f3* b ) + 7] + < - f |cu 2 w 3 ) 



(4.26) 



tt(£ + 1) 2(£+l) 



5 2 Ue(a a - /3£) + 



7T 



C(K» 



p / ii(a ai -a J2 ) I 1_ 

S 2 {2iea a + 27re + t]\u 2 uj 3 ) TT TT V 7r(€+1) 5+1 



,,1V ^2(2iea a + 7T^ + 27re + r7|cj 2 ^3) ^ 4 rf ^rs) ■ i N 

tt TT S 2 (i€(a ai + a a2 ) + 27re + r/la^a^^Qietaa! - a a2 ) + 7r|u;2^3) 

ai<a 2 e =± ^(^(aai + «a 2 ) + ^ + 2vre + T) \lv 2 lv 3 ) S 2 {ie{a ai - a a2 ) + vr(£ + l)|a; 2 a;3) 
tt tt r 2 (iea a + i/i £ - | + 7re|u;2^3)r2(-^a a + 77 - f - ?re + iji e \uj 2 uj 3 ) 
11 11 r 2 (iea a - i/i e + 7r£ + | + 7re|u;2^3)r2(-«ea a + ?7 + 7r£ + |-7re- i/i^c^s) 



x 

a<=A e=± 



rrrr ; i 2 yz,tta a -t zvie -r vi \m 2 , z,lu 3 )l 2 yz,tta a -+- zvie -|- vi ^ -|- i;|^ 2 , z,m 3 ) , , 

11 11 \l V J 0.4 en. 4- 9.Tfl/,) n 9y,)„^F„/'9iVn, -4- -4- <7r£k)„ 9i,iJ ' 

aGA e=± 

iL=±V r 2 (2* e a a + 2vre + 2^K 2u 3 )T 2 (2iea a + 27re + 2r? + 7r£|cu 2 , 2cu 3 ) ' 
The magnetization on a site m is given by 

({0}|{0}) 

(O = • (4-28) 

({0}|{0}) 

iVote. In paper j^J the authors considered the following vacuum expectation value for 
finite XXZ chain with double boundaries at massive regime, 



{vac 


e^(Cr)---$e^(a)e F *e G $ e]V 


(Ov)---$ ei (Ci)e G * 


vac) 


{vac 


e F e F *e G e G * 


vac) 



which is free from a difficulty of divergence. However it's physical meaning is not clear. 
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A Vertex Operators 

Here we summarize the bosonizations of the vertex operators || . 
Let us set free bosons b(t)(t G K) which satisfy 

sh I — I sh(7rt)sh 

W), Kt')\ = K2 { H , t fi 2 W + 1'). (a.i; 



Let us set a(t) by 



6( ( )sh^il = a (i )sh^. (A.2) 



The bosonization of the type-I vertex operators is given by 

$ + (/5) = £/(/?), (A.3) 
$_(/?) = /" da : U(P)U(a) : 

r l^ri) + WTT)J r l ^ny + 2^Ti)J' (A - 4) 



where we have set 



Kt) Jat ^\ . r - Tf . .( f 00 I'd) , 



U(a) =: exp ( - / ^dt :, U(a) =: exp f / ^e^dt ) : . (A.5) 



The bosonization of the type-II vertex operators is given by 

= 7(0), (A.6) 
= J da: V((3)V(a) : 

x r ?i^-lVf-^^ + --V (A.7) 
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where we have set 

= : (.L ss^J : - v(a) ; l" / srp 6 *"*; : • < A - 

Here the integration contours are chosen as follows. The contour Cj is (—00,00). The 
poles 

«-/?=— + mr(£ + l)i, (n £ N) (A.9) 

Z 

°f r + 2(|+T)) are a bove (7/ and the poles 

a-/3 = -— -rwr(f + l)i, (n G N) (A.10) 

Z 

°f r ^— ^^fj + 2(|+i) ) are below C/. The contour C/j is (—00, 00) except that the poles 

a-P= -— + nir&, (n e N) (A.ll) 

Z 



of T — ^ j are above C/j and the poles 

a-/3= rwrfi, (neN) (A. 12) 

Z 

of r (-*^ - ^) are below C/j. 

B Multi Gamma functions 

Here we summarize the multiple gamma and the multiple sine functions. 
Let us set the functions Y\{x\jS)^ r 2 (x\uii, u^) and T 3 (x\uji, 002, w 3 ) by 



logr 1 (x\u)+>yB u (x\u) = / — -e-** - &K ' (B.l) 



e -xt log(-*) 
c 27rit 1 - e~ ut 



logr 3 (x|wi, w 2 , ^3) + ■^B 33 (x\uj 1 , u 2 , u 3 ) = J 



dt _ e _ xt log(-t) 



c 2mt (1 - e-^*)(l - e-<" 2t )(l - e^*) ' 

(B.3) 

where the functions Bjj(x) are the multiple Bernoulli polynomials defined by 

±r xt 00 

— = V] — B rn (a;|wi • • -w r ), (B.4) 
1 ^—^ n\ 



m=i(^«-i) 



l J = iv ' n=0 
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more explicitly 

Bu(x\u) = (B.5) 
B 22{xW) = J^-(± + lWl + l(^ + ^Y (B.6) 

Here 7 is Euler's constant, 7 = ]im n _> 00 (l + | + | + l~ ^ — logn). 

Here the contor of integral is given by 




Contour C 



Let us set 



We have 



S 1 (x\u;) = — (B.7) 

Ii(u; — x|c<j)ri(a;|c<j) 

c / 1 \ r 2(^i + ^2 -x\u}i,u} 2 ) fu fi x 

5 2 (x wi,w 2 ) = w~/ i x , B.8 

r 2 {x\u 1 ,U2) 

5 3 (x|wi,w 2 ,W3) = Tw — 1 ^fth V ( R9 ) 

r 3 (o;i + lo 2 + uj 3 - x\u 1 ,u 2 , a; 3 )r 3 (x|a;i, u 2 , uj 3 ) 



T 1 (x\u) = e ^-^ loga, ^M , S^w) = 2sin(7nr/cu), (B.10) 

V 27T 



r 2 (a; + u; 2 ) _ 1 S , 2 (x + a; 2 ) _ 1 Ti(a; + cj|cj) 

r 2 (x|a;i,a; 2 ) Ti(a;|w 2 )' S , 2 (x|a;i, a; 2 ) S'i(a;|a; 2 )' r^xja;) 



= x.(B.ll) 



r 3 (x + u;i|a;i,^2,^3) _ 1 S 3 (x + u^Ui,^,^) _ 1 

r 3 (x|a;i,a; 2 ,a; 3 ) r 2 (x|a; 2 , a; 3 ) ' S , 3 (x|a;i, a; 2 , a; 3 ) S , 2 (x|a; 2 , a; 3 ) ' 
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r shfx — ^ 1+ ^ 2 )f dt 



27T 

S 2 (x\lu 1 lu 2 ) = ^=x + 0(x 2 ), (x->0). (B.14) 



7TX 7TX 

S 2 (x\uiUJ 2 )S 2 (— x\l)ilj 2 ) = — 4sin — sin — . (B.15) 

LOi L0 2 



21 



